Two Meson Systems with Ginsparg- Wilson Valence Quarks 
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Unphysical effects associated with finite lattice spacing and partial quenching may lead to the 
presence of unphysical terms in chiral extrapolation formulae. These unphysical terms must then be 
removed during data analysis before physical predictions can be made. In this work, we show that 
through next-to-leading order, there are no unphysical counterterms in the extrapolation formulae, 
expressed in lattice-physical parameters, for meson scattering lengths in theories with Ginsparg- 
Wilson valence quarks. Our work applies to most sea quark discretization, provided that chiral 
perturbation theory is a valid approximation. We demonstrate our results with explicit computations 
and show that, in favorable circumstances, the extrapolation formulae do not depend on the unknown 
constant Oyiix appearing at lowest order in the mixed action chiral Lagrangian. We show that the 
1 = 1 KK scattering length does not depend on CMix in contrast to the 1 = 3/2 Kir scattering 
length. In addition, we show that these observables combined with /k//tt and the 1 = 2 irn 
scattering length share only two linearly independent sets of counterterms, providing a means to 
test the mixed action theory at one lattice spacing. We therefore make a prediction for the 1 = 1 KK 
scattering length. 

PACS numbers: 12.38.Gc 



I. INTRODUCTION 

There is currently a tension in lattice simulations of 
QCD phenomena between the need for quarks obeying 
chiral symmetry on the lattice, and the need for quark 
masses light enough that one is in the chiral regime. 
This tension occurs because quark discretization schemes 
which obey chiral symmetry on the lattice, such as do- 
main wall fermions [l|, 0, H[ or overlap fermions 0, V 
both of which satisfy the Ginsparg- Wilson relation 
are numerically expensive to simulate. On the other 
hand, Wilson fermions [§] or staggered fermions [l(J, [ll[ 
are faster but violate chiral symmetry at non-zero lattice 
spacing. 

One way of resolving this tension is to recognize that 
the most computationally intensive stage of a fully dy- 
namical simulation is the evaluation of the quark de- 
terminant. This determinant is associated with the sea 
quarks and is a component of the probability measure 
on the space of gauge field configurations. This observa- 
tion has l ong been the motivation for partial quenching 
(PQ) [H, [U — sea quark masses are taken to be larger 
than valence quark masses so that the sea quarks are 
more localized and the determinant is easier to compute. 
The notion of a "mixed action" (MA) simulation takes 
this line of reasoning one step farther [3, [ll| . A mixed 
action simulation uses different quark discretizations in 
the sea and valence sectors. In this case, the valence 



quarks can be chosen to obey the Ginsparg- Wilson rela- 
tion so that they enjoy chiral symmetry at finite lattice 
spacing. The numerically expensive sea quarks, on the 
other hand, can be chosen to be inexpensive Wilson or 
(rooted) staggered sea quarks, for example. 

There have recently been a significant number of mixed 
action lattice QCD simulations, [H, [U GJ, GJ, E EH, 
HI, [H, H3, HH, HH, the majority of which have employed 
domain wall valence fermions on the publicly available 
MILC lattices [27j . 1 The effective theories appropri- 
ate for mixed action simulations were originally devel- 
oped in Refs. 0, [H| for Ginsparg- Wilson (GW) valence 
fermions on Wilson sea fermions, and later for GW va- 
lence fermions on staggered sea fermions |3(J; these the- 
ories have received considerable theoretical attention re- 
cently [U, H3, [H, [H, [H, H6j in response to the numerical 
interest. 

In this article, we study aspects of the chiral per- 
turbation theories appropriate for mesonic processes in 
mixed action simulations with Ginsparg- Wilson valence 
quarks. The chiral properties of the Ginsparg- Wilson 
valence quarks are central to our work, so we adopt the 
convention that when we refer to a mixed action simu- 
lation, we imply that the valence quarks are Ginsparg- 
Wilson. We work consistently at next-to-leading order 
(NLO) in the effective field theory expansion, which is 
a dual expansion in powers of the quark mass, m q and 
the lattice spacing, a. At this order, one can view cur- 
rent lattice simulations as being methods of computing 
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the values of certain coefficients in the NLO chiral La- 
grangian, which for mcsonic quantities is known as the 
Gasser-Leutwyler Lagrangian [33, [38[. This is because 
lattice simulations are performed at quark masses larger 
than the physical quark mass, so that lattice data must 
be fit to formulae computed in chiral perturbation theory. 
These fits determine the unknown coefficients occurring 
in the chiral formulae, known as low energy constants 
(LEC)s, so that the chiral expression can then be used 
at the physical values of the meson masses and decay 
constants to predict the results of physical experiments. 
Frequently there are non-physical operators in the NLO 
chiral Lagrangians describing discretized fermions. For 
example, there are of order 100 operators in the NLO 
staggered chiral Lagrangian [39| compared to order 10 
in the Gasser-Leutwyler Lagrangian. These unphysical 
operators lead to unphysical terms in chiral extrapola- 
tion formulae which must somehow be removed to make 
physical predictions. One might think that this will also 
be an issue in mixed action chiral perturbation theory 
(MAxPT), since there are certainly many additional op- 
erators at NLO. However, we show that mixed action 
lattice simulations of mesonic scattering lengths do not 
depend on any unphysical operators at NLO, if these 
scattering lengths are expressed in terms of the pion mass 
measured on the lattice and the decay constant measured 
on the lattice [40(. For linguistic brevity, we will refer 
to the pion mass measured on the lattice as the lattice- 
physical pion mass and similarly for the decay constant. 

Each choice of sea quark discretization leads to a differ- 
ent MAxPT. For example, tree level shifts of the masses 
of mesons composed of two sea quarks, or of one valence 
and one sea quark, are different for staggered sea quarks 
and for Wilson sea quarks. Therefore, one may think that 
the chiral extrapolation formulae depend on the nature 
of the sea quark discretization. We show that, at NLO, 
the only difference between the extrapolation formulae 
is in the leading order mass shifts of the mesons com- 
posed of two sea quarks. Thus, once these mass shifts 
are known, one can use the same extrapolation formulae 
for different sea quark discretizations. 2 In fact, any sea 
quark discretization will do provided, 3 firstly, that QCD 
is recovered in the continuum limit, and secondly, that 
the sea-sea mesons may be described at leading order 
by chiral perturbation theory with the usual kinetic and 
mass terms, or that the non-locality of the appropriate 
chiral perturbation theory is correctly captured by the 
replica method 42]. 4 In addition, we have assumed that 



2 For staggered fermions, these mass corrections to the mesons are 
well know [4l| |. However, these effects are less well determined 
for Wilson fermions. 

3 In this paper, we restrict ourselves to isospin symmetric masses 
in the valence and sea sectors. 

4 We have in mind the current discussion regarding whether rooted 
staggered fermions become QCD in the continuum limit. There 
is growing numerical and formal evidence lending support to the 



the quarks of the sea sector are only distinguished by 
their masses, so that, for example, the same discretiza- 
tion scheme has been used for all sea quarks, and we 
assume that chiral perturbation theory itself is a valid 
approximation. 

There are still various challenges facing mixed action 
simulations. Mixed action simulations always violate uni- 
tarity at finite lattice spacing. 5 In MA%PT and PQxPT, 
the most severe unitarity violations are encoded in hair- 
pin propagators of flavor neutral mesons. We point out a 
simple parametrization which allows a convenient book- 
keeping of these unitarity violating effects. Addition- 
ally, the value of the new constant, Civiix, that appears 
in the LO mixed action Lagrangian, is currently un- 
known. This term leads to an additive lattice spacing 
dependent mass shift of "mixed" mesons consisting of 
one valence and one sea quark. This causes a mismatch 
of the meson masses composed of different quarks but 
does not play a role in the well known enhanced chiral 
logarithms pjj 0, 5|| or the enhanced power-law vol- 
ume dependence of two-hadron states . 6 In addition, 
the value of this constant is presumably different for each 
sea quark discretization. However, we show that under 
favorable circumstances physical quantities such as scat- 
tering lengths do not depend on this constant, as was 
first noticed in Ref. [341 ]. 

Finally, to demonstrate our arguments, we determine 
various NLO formulae for use in chiral extrapolation of 
certain mcsonic quantities. We have computed the KK 
and Kit scattering lengths in SU(6\3) MAxPT. In ad- 
dition, for completeness, we present the irir scattering 
length in SU(4\2) and SU(6\3) MAxPT, which were 
computed in Ref. 34] as well as the ir and K meson 
masses and decay constants, which were first computed 
in Refs. 0, EH [3(3] , but we express these quantities in 
terms of the PQ parameters we introduce in Sec. Ill Al We 
conclude with a discussion of our results, and in particu- 
lar show that among the three meson scattering lengths 
mentioned above and the quantity fx/ fir there are only 
two linearly independent counterterms at NLO, which are 
the corresponding physical counterterms of xPT. There- 



hope that rooted staggered fermions are in the same universal- 
ity class as QCD. We refer the reader to Ref. for current 
summary of the issues. 

5 For arbitrarily small lattice spacings, the differences arising from 
the different lattice actions will become negligible, which prac- 
tically means smaller than the statistical and systematic un- 
certainties for a given observable, and these unitarity violating 
terms will no longer be important (assuming the quark masses 
are tuned equal). This of course, also implies that a MA effective 
field theory description will no longer be necessary, however for 
MA lattice simulations today and the foreseeable future, MAxPT 
is the necessary tool for controlling extrapolations to the physical 
point. 

6 In this article we do not discuss the observed negative norm issues 
involving scalar meson correlators [4ll . [4TJ , but these can at least 
be qualitatively understood with the appropriate effective field 
theory methods [illil. 
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fore, these four processes provide a means to test the MA 
formalism with only one lattice spacing. In the appen- 
dices we collect the various formulae which are necessary 
for the chiral extrapolations of the quantities we discuss 
in this paper. 



II. MIXED ACTION EFFECTIVE FIELD 
THEORY 

We will not give a thorough introduction to mixed ac- 
tion or partially quenched theories here. We will simply 
give a brief review to introduce our notation and power 
counting. For a good introduction to MA theories we 
refer the reader to Refs. [3, HH, [3(3], and for PQ theories 
to Refs. [HiH. 



A. Mixed Actions at Lowest Order 

To construct the appropriate Lagrangian to a given 
order, one must specify a power counting. As mentioned 
above, X-PT is a systematic expansion about the zero 
momentum, zero quark mass limit, for which the small 
expansion parameter is 
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where m 2 cx m q . For effective theories extended to in- 
clude lattice spacing artifacts, one must include an addi- 
tional small parameter. 7 We will be interested in theories 
for which the leading sea quark lattice spacing depen- 
dence is 0(a 2 ), such as staggered, 0(a)-improved Wil- 
son (55j . twisted-mass at maximal twist [56], or chiral 
fermions. Therefore, we shall denote the small parame- 
ter counting lattice spacing artifacts to be 



n 2 A 2 

a 1V QCD 



(2) 



and we shall work consistently in the dual expansion to 

0(4) , 0(44) , o(si). (3) 

At leading order (LO) in the quark mass expansion, the 
mixed action Lagrangian is simply given by the partially 
quenched Lagrangian 

C = ^str (d^d^) + ^pstr (m 9 E t + Em*) , (4) 



where we use the normalization / ~ 132 MeV, and 
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The matrices M and M contain bosonic mesons while \ 
and x contain fermionic mesons with one ghost quark 
or antiquark. To be specific, we will discuss the theory 
with 3 valence (and ghost) quarks and 3 sea quarks, for 
which 8 
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The upper N v x N v block of M contains the usual mesons 
composed of a valence quark and anti-quark. The lower 
N s x N s block of M contains the sea quark-antiquark 
mesons and the off-diagonal block elements of M contain 
bosonic mesons of mixed valence-sea type. 

For MA theories there are two types of operators we 
need to consider at LO in e 2 . There are those which 
modify the sea-sea sector meson potential, which we shall 
denote as U sea , and those which modify the mixed meson 
potential, which we shall denote as Uys-, such that the 
Lagrangian, Eq. Q is modified by the additional terms 
(following the sign conventions of Ref. (30j). 



I->MA = ~ a 2 [Usea — Uy S 



(7) 



We shall not specify the form of the sea-sea meson poten- 
tial, U sea , but only note again that at the order we are 
concerned with, we only need to know how the masses 
of the sea-sea mesons are modified at LO in s 2 which 
have been discussed, for example, in Refs. [3, Ho|. The 
other important thing to know is that the structure of 
Uys is independent of the type of sea quark and is given 



by 



s is ma 

MM 



Uy S = C M iK str(r 3 S T 3 , 



(8) 



where the flavor matrix T3 is a difference in projectors 
onto the valence and sea sectors of the theory, 



T 3 = V S - Vy = diag(-/y,7 s , -I v ) 



(9) 



This operator leads to an additive shift of the valence- 
sea meson masses, such that all the pseudo-Goldstonc 



7 The general procedure p$2j is to construct the continuum 
Symanzik quark level effective theory for a given lattice ac- 
tion [53L and then build the low energy effective theory with 
spurion analysis on this continuum lattice action. 



8 For staggered sea quarks, each sea quark label implicitly includes 
a taste label as well. For example, <j> u j is a 1 X 4 vector in taste- 
space. 
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mesons composed of either valence quarks, v, sea quarks, 
s, or both have LO masses given by 9 



TO « lt , 2 = B (m Vl + m V2 ) , 
rh 2 vs = B (m v +m s ) + a 2 A M ix : 
B (m Sl + m S2 ) + a 2 A sec 



with A s 



determined by U sea and 10 

16 C*Mix 
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The most severe and well known unitarity violating 
feature of PQ and MA theories is the presence of double 
pole propagators in the flavor neutral mesons (T^ |. In 
particular the momentum space propagators between two 
mesons composed of valence quarks of flavors a and b 
respectively are given by 11 



theories, when these quantities are zero, the theory re- 
duces to an unquenched theory. For MA theories, when 
one tunes these parameters to zero, one tunes the double 
pole structure of the flavor-neutral meson propagators to 
zero up to higher order corrections, and thus has the most 
QCD-like scenario for a MA theory (we note that there 
is still a mismatch in the mass of the mixed meson, com- 
posed of one valence and one sea quark, Eq. (|10p . from 
the others) . We therefore introduce the partial quenching 
parameters, 12 



X 2 _ - 2 _ 2 
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a 2 A, 



a 2 A sea 



(14) 



where the dots denote higher order corrections to the 
meson masses. We will now move on to discuss the gen- 
eral structure of mixed action theories for arbitrary sea 
quarks at the next order. 
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where rhx is the mass of the 7y sea -field, 



B. Mixed Action xPT at NLO 

It was shown in Ref. [HI that the I = 2 tttt scattering 
length at NLO, expressed in terms of the bare parameters 
of the chiral Lagrangian, is 13 



m 2 x 



1™% 



(13) 



When the valence quark masses are equal, either in the 
isospin limit of light quarks or for the same flavor, a = b, 
the above propagator acquires a double pole. It is these 
double poles which lead to the well-known sicknesses of 
PQxPT, such as the enhanced chiral logs [H, 0, [45[ , and 
enhanced power-law volume dependence of two-particle 
states 46]. Here, we introduce what we call "partial 
quenching parameters," which are a difference in the 
quark masses for PQ theories and, more generally for 
MA theories, a difference in the masses of mesons com- 
posed of two sea quarks and two valence quarks. For PQ 



Here and throughout this article, we use tildes over the masses to 
indicate additive lattice spacing corrections to the meson masses. 
For a twisted mass sea 1591 . one must keep separate track of 
the neutral and charged mesons as they receive a relative O(o?) 
splitting [57l |. similar to the various taste mesons for staggered 
fcrmions. 

In mixed action theories, there are additional hairpin interactions 
proportional to the lattice spacing which arise from unphysical 
operators in the theory, similar to the lattice spacing dependent 
hairpin interactions in staggered xPT 58]. For O(a) improved 
Wilson fermions and staggered fermions, these effects are higher 
order than we are concerned with in this paper l.'i'J , which is also 
true for twisted mass fermions at maximal twist. For Wilson 
and twisted mass fermions (away from maximal twist), these 
effects appear at the order we are working, and must be included. 
We assume for the rest of the paper that the sea quark scaling 
violations are O(o?) or higher. 
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Let us point out some features of this expression which 
are relevant from the point of view of chiral extrapo- 
lations. Equation (fl~5| depends on the mass rhj u of a 
mixed valence-sea meson and consequently the expres- 
sion depends on the value of the parameter Cmix- In 
addition, there is a dependence on the unphysical un- 
knowns t'pQ (n) and £' a 2 (/i) , as well as the decay constant 
and chiral condensate in the chiral limit, / and Bq- Thus, 



For a staggered sea it is the taste-identity meson masses which 
enter these PQ parameters |30H (which have been measured [4l|) 
while for a twisted mass sea, it is the neutral pion mass. 
13 Here, we show the scattering length for a two-sea flavor theory. 
In the appendix, we also list the result for the three-sea flavor 
theory. However, the following discussion of the counterterm 
structure of the NLO Lagrangian is independent of the number 
of sea flavors. 
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Eq. (|15[) depends on three unphysical unknown param- 
eters and three physical unknown parameters, £' vn (fj,),f 
and Bq. One must fit all unknown parameters to extrap- 
olate lattice data but only three are of intrinsic interest. 

In terms of lattice-physical parameters, the same scat- 
tering length becomes 



,1=2 
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Notice that this expression does not depend on the mixed 
valence-sea mesons, and, in fact, the only unknown terms 
in the expression is the physical parameter £ n „(fj,), and 
the sea-sea meson mass shift in A| M , Eq. (fTX)]) . which is 
already determined for staggered sea-quarks. Thus, chi- 
ral extrapolations using the formula Eq. (|16p requires 
fitting only one parameter (two for non-staggered sea 
quarks), in contrast to chiral extrapolations using the 
scattering length expressed in terms of the bare parame- 
ters, Eq. (fT5|) . Our goal in this section is to understand 
the origin of this simplification, and under what circum- 
stances we may expect similar simplifications to occur in 
other processes. To do so, we must discuss the structure 
of the NLO terms of the MA^PT Lagrangian. 

The symmetry structure of the underlying mixed ac- 
tion form of QCD determines the NLO operators in the 
mixed action chiral Lagrangian through a spurion anal- 
ysis. However, the symmetries enjoyed by the valence 
quarks are different to the symmetries of the sea quarks 
in a mixed action theory. In particular, we only consider 
GW valence quarks which have a chiral symmetry; the 
numerically cheaper sea quarks typically violate chiral 
symmetry. Thus, it is helpful to consider spurions aris- 
ing from the valence sector separately to the spurions of 
the sea sector. 

The valence sector only violates chiral symmetry ex- 
plicitly through the quark mass. Therefore, at NLO, 14 
the purely valence spurions are identical to the spuri- 
ons in continuum, unquenched chiral perturbation the- 
ory, and so the valence- valence sector of the NLO mixed 
action chiral Lagrangian is the Gasser-Leutwyler La- 
grangian. The sea sector is different. At finite lattice 
spacing, the sea sector has enhanced sources of chiral 
symmetry violation — for example, there are additional 
spurions associated with taste violation if the sea quarks 
are staggered, or in the case of a Wilson sea, the Wilson 
term violates chiral symmetry. Consequently, there are 



14 Lattice artifacts such as Lorentz symmetry violation lead to the 
presence of unphysical operators in the chiral Lagrangian. These 
operators will not be important in the following, as they are 
higher order in the chiral expansion for mesons |l5j (however 
they are relevant at for baryons [3ll0 . 



additional spurions in the sea sector. Of course, these 
spurions must involve the sea quarks and must vanish 
when the sea quark fields vanish. 15 Nevertheless, scat- 
tering amplitudes expressed in terms of lattice-physical 
parameters do not explicitly depend on the lattice spac- 
ing a, as we will now discuss. 

In this paper, we work consistently to NLO in the 
MA%PT power counting which we have defined in 
Eq. ([3]). At this order, the NLO operators in the La- 
grangian are only used as counterterms; that is, at NLO 
one only computes at tree level with the NLO opera- 
tors. Since the in/out states used in lattice simulations 
involve purely valence quarks, we must project the NLO 
operators onto the purely valence quark sector of the the- 
ory. Consequently, all of the spurions which involve the 
sea quark fields vanish. Since the remaining spurions in- 
volve the valence quarks alone, we only encounter the 
symmetry structure of the valence quarks as far as the 
NLO operators are concerned. These spurions only de- 
pend on quark masses and the quark condensate itself, 
and so there can be no dependence on lattice discretiza- 
tion effects arising in this way. The exception to this 
argument arises in the case of double trace operators in 
the NLO chiral Lagrangian; in these cases the valence 
and sea sectors interact in a flavor-disconnected manner, 
unlike the operator in Eq. ([5]). If one trace involves a 
valence-valence spurion while the other involves a sea- 
sea spurion, then the trace over the sea may still con- 
tribute to a physical quantity, for example the meson 
masses and decay constants, see App. lATICl for explicit ex- 
amples. Note that the valence-valence operators which 
occur in these double trace operators must be propor- 
tional to one of the two operators present in the LO chi- 
ral Lagrangian, Eq. Thus for meson scattering pro- 
cesses, the dependence upon the sea quarks from these 
double trace operators can only involve a renormaliza- 
tion of the leading order quantities / and Bq. Both the 
explicit sea quark mass dependence and the explicit lat- 
tice spacing dependence are removed from the scattering 
processes expressed in terms of the lattice-physical pa- 
rameters since they are eliminated in favor of the decay 
constants and meson masses which can simply be mea- 
sured on the lattice. We therefore conclude that when 
expressed in lattice-physical parameters, there can be no 
dependence upon the sea quark masses leading to un- 



15 Not all the lattice spacing dependence may be captured with spu- 
rion analysis. There are O(o?) operators at the quark level which 
do not break chiral symmetry, for example, 0( 6 ' = a 2 Q I/P Q. 
This operator leads to an a 2 renormalization of all the low energy 
constants (LEC)s of the low energy theory. Because this opera- 
tor does not break any of the continuum QCD symmetries, it can 
not be distinguished through spurion analysis |l5|. Of course, an 
operator of this form is present when the QCD Lagrangian is run 
from a high scale (say, the weak scale) down to the scale of the 
lattice, so its effects could in principle by accounted for by per- 
forming a perturbative matching computation between the QCD 
effective Lagrangian at the scale fi = a -1 and the lattice action. 
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physical PQ counterterms and similarly there can be no 
dependence upon an unphysical lattice-spacing countert- 
erm. 

Let us present another more physically intuitive argu- 
ment concerning the absence of sea quark mass depen- 
dence in meson scattering processes. To do so, we must 
digress briefly on tttt scattering in SU(3) chiral perturba- 
tion theory. The strange quark mass m s is a parameter 
of SU(3) xPT, and so one would expect that the tttt 
scattering length includes analytic terms involving m s . 
However, we can consider a theory in which the strange 
quark is heavy, so that we may integrate it out; we must 
then recover SU{2) xPT. Chiral symmetry forces any m a 
dependence in the analytic terms of the tttt amplitude to 
occur in the form m^m 2 K . But the only counterterm in 
the on-shell SU{2) scattering amplitude (Eq. (fl6j) with 
A JU = 0) is proportional to m%. It is not possible to ab- 
sorb m^rn 2 K into so there can be no m s dependence 
in the SU(3) tttt scattering amplitude. This is indeed the 
case, as was observed in Ref. [59(. 

Now, let us return to the PQ and MA theories. For 
the purposes of this discussion, we can ignore the flavor- 
neutral and ghost sectors, reducing our theory from an 
SU(6\3) theory to an SU(6) theory. The sea quark de- 
pendence of this 5(7(6) chiral perturbation theory is anal- 
ogous to the m s dependence of SU(3) x?T in our tttt 
scattering example (as in this process the strange quark 
of SU(3) only participates as a sea quark). A similar de- 
coupling argument tells us that the sea quark masses can- 
not affect processes involving the valence sector provided 
one uses the analogues of on-shell parameters which are 
the lattice-physical parameters. We conclude that there 
can be no analytic dependence on the sea quark masses 
in a mesonic scattering amplitude. Further, these ar- 
guments only depend upon the chiral symmetry of the 
valence quarks and thus also apply to the lattice spacing 
dependence. 

Now, we shall make these arguments concrete by ex- 
plicit computations. The NLO Lagrangian describing the 
valence and sea quark mass dependence is the Gasser- 
Leutwyler Lagrangian with traces replaced by super- 
traces: 

£-gl = L\ [str (c^Ec^Et)] 2 

+ L 2 str (c^Ec^E*) str (d^d v ^) 

+ L 3 str (d^d^d^Ed^) 

+ 2B L 4 str (c^Ec^E*) str (rriqY^ + EmJ) 

+ 2P L 5 str [c^Ec^E* (m g E t + Em*)] 

+ AB 2 L 6 [str(m 9 E t + Emt)] 2 

+ 4:B 2 L 7 [str (mjE-Etmg)] 2 

+ 4B 2 L s str (m g E t m 9 E t + Em* £ m t) . (17) 



Having a concrete expression for the Lagrangian, 16 we 
can easily show explicitly how the sea quark mass de- 
pendence disappears. The key is that when constructing 
NLO correlation functions of purely valence quarks, we 
can replace the mesonic matrix $ in the NLO Lagrangian 
by a projected matrix 

$ -> P V <S>P V (18) 

where Py is the projector onto the valence subspace. 
Therefore the matrix E has an expansion of the form 

E = 1 + Py^Pv + ■■■ (19) 

Now, insert this expression into Eq. (jTTJ) , and consider 
only the terms involving non-zero powers of In the 
single trace operators, the projectors remove any depen- 
dence on the sea quark masses. There is still sea quark 
mass dependence remaining in the double trace operators 
proportional to L4 and L$ given by 

5C GL = 4P U str (c^EPy d^E+Py) str(m 8 ) 

+ 16P^ L 6 str (m q ^P v + PyEmJ) str(m g ). (20) 

However, these operators simply shift / and B$ 

f 2 -> f + 321-4 B str(m g ) (21) 
f 2 B Q -» f 2 B Q + 64L 6 Bl str(m 9 ). (22) 

Since the parameters / and Bo are eliminated in lattice- 
physical parameters in favor of the measured decay con- 
stants and meson masses, we can remove the dependence 
of scattering lengths on the sea quark masses by working 
in lattice-physical parameters. In an analogous way, we 
can remove all the explicit lattice spacing dependence. 
The general MA Lagrangian involving valence-valence 
external states at 0(e 2 n e 2 ) can be reduced to the fol- 
lowing form 

SCma = a 2 L d a2 str (c^EPy d^P v ) 

x str (f(PsZPs)f'(Ps^Ps)) 
+ a 2 L™ q str (m q P v ^P v + PyEPymJ) 
x str (. 9 (P S EP S ) g'(Ps^Ps)) + h.c, (23) 



The generators of the PQ and MA theories form graded groups 
and therefore lack the Cayley-Hamilton identities of SU (N) the- 
ories. Therefore, PQ and MA theories have additional oper- 
ators compared to their counterparts. For example, the 
C(p 4 ) Lagrangian has one additional operator as compared to 
the Gasser-Leutwyler Lagrangian. However, we do not need to 
consider the effects of this operator in our analysis as it has been 
shown that it can be constructed such that it does not contribute 
to valence quantities until 0(p 6 ) ,60). This is not generally the 
case, as is demonstrated by various examples in the baryon sec- 
tor [m, m, m, m, m, m, [S3, m, m, n^i . 
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where the /'s and g's are functions dependent upon the 
sea-quark lattice action. These then lead to renormaliza- 
tions of the LO constants, 

f -> f + 8a 2 L% str (/(P S SP S ) f'(P s tfP s )) , 

f 2 B -► / 2 B + 4a 2 i™' str ( ff (P s EP s ) 5 '(P s £tp s )) , 

(24) 

and just as with the sea quark mass dependence, express- 
ing physical quantities in terms of the lattice-physical 
parameters removes any explicit dependence upon the 
lattice spacing in mesonic scattering processes. 

Together, these results show that at NLO, the only 
counterterms entering into the extrapolation formulae 
for mesonic scattering lengths are the same as the coun- 
terterms entering into the physical scattering length at 
NLO. This lack of unphysical counterterms is desirable 
from the point of view of chiral extrapolations, but it also 
has another consequence. Loop graphs in quantum field 
theories are frequently divergent; there must be a coun- 
terterm to absorb these divergences in a consistent field 
theory. Since there is no counterterm proportional to 
a 2 or the sea quark masses, loop graphs involving these 
quantities are constrained so that they have no diver- 
gence proportional to a 2 or the sea quark masses. This 
further reduces the possible sources of sea quark or lat- 
tice spacing dependence. For example, mixed valence-sea 
meson masses have lattice spacing shifts, so there can be 
no divergence involving the valence-sea meson masses. In 
some cases this constraint is strong enough to force the 
entire valence-sea mass dependence to cancel from scat- 
tering lengths expressed in lattice-physical parameters. 
If this occurs, then the scattering length will not depend 
on the unknown constant CmIx- 



Dependence upon sea quarks 

Now, let us move on to discuss how the NLO extrap- 
olation formulae depend on the particular sea quark dis- 
cretization in use. At NLO in the effective field theory 
expansion, mesons composed of one or two sea quarks 
only arise in loop graphs. In particular, the valence-sea 
mesons can propagate between vertices where they inter- 
act with valence-valence mesons; these interactions in- 
volve the LO chiral Lagrangian augmented with the mix- 
ing term a 2 Uvs- Because the mixing term is universal, 
these interaction vertices are the same for all discretiza- 
tion schemes provided LO chiral perturbation theory is 
applicable. The sea-sea mesons only arise at NLO in hair- 
pins; therefore, they are only sensitive to the quadratic 
part of the appropriate LO chiral Lagrangian on the sea- 
sea sector. Thus, we see that our NLO extrapolation 
formulae only depend on the LO chiral Lagrangian to 
quadratic order in the sea-sea sector and the LO chi- 
ral Lagrangian (with the mixing term) in the valence-sea 
sector. Together, we see that the condition we require 
on the sea quark discretization is that the sea-sea sector 



alone should be described by chiral perturbation theory 
at LO, and that the constant C*Mix should not be so large 
that its explicit violation of chiral symmetry overwhelms 
the dynamical violation of chiral symmetry. Non-locality 
which is described by the replica trick does not present 
a problem since at the level of perturbation theory the 
analytic continuation required by the replica method is 
trivial. 

Note that the impact of using different sea quark 
discretizations in our work is only at the level of the 
quadratic Lagrangian. Therefore, the same NLO extrap- 
olation formulae can be used to describe simulations with 
different sea quark discretizations, provided that the ap- 
propriate mass shifts are taken into account. In the case 
of staggered sea quarks, the sea-sea mass splitting which 
occurs in the MA formulae is that of the taste-identity, 
which has been measured [4l[, and for the coarse MILC 
lattices, is given by 

a 2 A seQ = a 2 A 7 ~ (450 MeV) 2 , (25) 

for a ~ 0.125 fm. These mass shifts can only ap- 
pear through the hairpin interactions at this order. 
These terms will generally be associated with unphysical 
MA/PQ effects which give rise to the enhanced chiral 
logarithms as well as additional finite analytic depen- 
dence upon the sea-sea as well as valence-valence meson 
masses (and their associated lattice spacing dependent 
mass corrections). The exception to this is the depen- 
dence upon the 77- mass. As can be seen in Eq. (|T2"j) . the 
only way the 77-mass dynamically enters processes involv- 
ing external pions and kaons through ©(e^e 2 ) is via the 
mass of the sea-sea rj, Eq. (p~3|) . The other way these 
discretization effects enter MA formulae is through the 
mixed valence-sea meson masses, Eq. (ITUl) . Currently, 
this mass shift, a 2 AMix, is not known for any type of sea 
quark discretization. This is one of the more important 
MA effects, because it enters many quantities of interest 
at the one-loop level, for both mesons and baryons, and 
thus to perform chiral extrapolations properly, this mass 
splitting must be taken into account. 

Mixed actions at NNLO 

It is important to note that these conclusions will not 
hold at NNLO in the effective field theory expansion. At 
this order, NNLO terms in the effective Lagrangian will 
introduce a 2 shifts of the Gasser-Leutwyler parameters 
themselves. In simulations which are precise enough to 
be sensitive to NNLO effects in chiral perturbation the- 
ory, these effects would have to be removed. In addition, 
there will be new effects which can not be absorbed into 
the Gasser-Leutwyler parameters, but are truly new lat- 
tice spacing artifacts. The simplest example to under- 
stand is to consider how the pion mass is modified at 
0(ef n s 2 ) in a MA theory with staggered sea quarks [70] |. 

Briefly, there will be contributions to the pion mass 
which break taste, arising for example from the Gasser- 
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Leutwyler operator in Eq. (fl~7|) with coefficient Lq. The 
taste breaking contributions to the pion mass arise when 
the valence pion is contracted with the meson fields in one 
of the super-traces while the other super-trace is taken 
over sea-sea mesons which form a loop at this order, 



a means to test the mixed action formalism with only one 
lattice spacing. 18 

A. JkIU 



Sml(NNLO) 



64m 



m ■ 



B (m Sl + m S2 ) 



F,t 



{A-Kff 



x m 2 



s±s 2 ,t 



111 



(26) 



where N s = 1/4 is the factor one inserts according to 
the replica method to account for the 4 t,l -root of the sea 
quark determinant and n t counts the weighting of the 
mesons of various taste propagating in the loop. The 
staggered meson mass of flavor F, and taste t, is given 
at LO by HHHHzIl 



-2 



111 



sis 2 ,t 



B (m si +m S2 ) + a 2 A(£ t ) 



(27) 



These taste-breaking effects are unphysical and their as- 
sociated /i-dependence can only be absorbed by the ap- 
propriate unphysical lattice spacing dependent operators 
arising in the mixed action Lagrangian. This is simply 
one of many possible examples of how the continuum-like 
behavior of mixed action theories will break down. 



III. APPLICATIONS 

In this section, we discuss applications of these re- 
sults to some specific quantities of physical interest. 
There have been a number of recent lattice computa- 
tions [H GJ, El M, SH, [H, S3, HI H utilizing the 
scheme first developed by the LHP collaboration [72|, [73| 
of employing domain wall valence quarks with the pub- 
licly available MILC configurations. In particular, the 
NPLQCD collaboration has computed the I — 2 tttt scat- 
tering length [19(, 17 fx/ fir [23j and determined both the 
7 = 3/2 and I = 1/2 Ktt scattering lengths through 
a direct determination of the I = 3/2 Ktt scattering 
length 25:]. As we will demonstrate by explicit computa- 
tion, the I — 1 KK scattering length, together with the 
above three systems share only two linearly independent 
sets of counterterms, which are the physical countert- 
erms of interest. Therefore, these four quantities provide 



The pion and kaon decay constants were computed 
in a mixed action theory with staggered sea quarks in 
Ref. [3(3] • In Appendix \C\ we include the general form 
of these results for arbitrary sea quarks to NLO, which 
we express in terms of the PQ parameters we introduced 
in Eq. (fT4|) . We use these formulae to estimate the er- 
ror arising from the finite lattice spacing in the recent 
determination of fx/ f-x in Ref. (23[, in which the contin- 
uum xPT form of this quantity was used to extrapolate 
the lattice data to the physical point. The MA functional 
form of this quantity depends upon the mixed valence-sea 
meson masses, and so we can not make a concrete pre- 
diction of the error made in this approximation, as the 
mixed meson mass depends upon Cmix, Eq. (|10p . which 
is currently unknown. 19 Consequently we form the ratio, 



fK 




fK 














MA 




QCD 




fK 






U 






QCD 



(28) 



and in Eq. (|C4[) , we provide the explicit formula for this 
quantity, with the mass tuning used in Ref. [23| (m qs = 



in,, 



A 2 



2 A/ ~ (450 MeV) 2 ). In Fig. [U 
we plot this ratio as a function of the mixed meson split- 
ting in the range -(600 MeV) 2 < a 2 A Mix < (800 MeV) 2 . 
We take the value of L5 (/j.) from Ref. [23j , as their various 
fitting procedures produced little variation in the extrac- 
tion of Ls(fj,). This provides us with an indirect means at 
estimating the error in the extrapolation of the quantity 
/if //tt- As can be seen from Fig. [TJ reasonable values of 
a 2 AMix can produce deviations in fx/ fn on the order of 
5%. These deviations are important enough to include 
in the fitting procedure (although still within the confi- 
dence levels in Ref. [HI El])' but not significant enough 
to determine a 2 AMix directly from the data in Ref. [23j . 
One can also determine the size of the hairpin contribu- 
tions alone by setting a 2 Ajyiix = 0, and, as can be seen 
in Fig. [TJ these effects are a fraction of a percent for all 
values of the pion mass. 



17 In addition to the lattice spacing modifications of the 1 = 2 tttt 
scattering length computed in Ref. [3411 . the exponential fi- 
nite volume corrections to this quantity were also computed in 
Ref. [74)1 . It was found that for the pion masses in use today, 
these effects were not significant, being on the order of 1%. It is 
expected that the exponential volume dependence in the other 
scattering processes will be similar to that of the two-pion sys- 
tem, as in all cases the pion is the lightest particle and will dom- 
inate the long range (finite volume) effects. 



In Ref. [34T ]. it was argued that to all orders in perturbation 
theory, the unitarity violating features of MA and PQ theories 
do not invalidate the known method of extracting infinite vol- 
ume scattering parameters from finite volume correlation func- 
tions 75, 76, 77], for all "maximally stretched" two- meson states, 
i.e. the 1 = 2 tttt, I = 3/2 Kir and 1 = 1 KK systems. 
In Ref. |35| . this quantity was recently estimated by comparing 
the MA form of the pion form-factor to a MA simulation [l7l |. 
Unfortunately, only one of the lattice data points was in the 
chiral regime, so a precise determination of this quantity was 
not possible. 
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Lattice Spacing corrections to fit/fji 



A(f ) 




FIG. 1: We plot the ratio, A(f K /fn) defined in Eq. J28]) 
as a function of the unknown mixed meson mass splitting, 
-(600 MeV) 2 < a 2 A M ix < (800 MeV) 2 . The observed de- 
viation from the continuum x?T formulae is on the order of 
5%, which is important, but not significant enough to directly 
determine this unknown mass splitting from the MA lattice 
data of fi</ fn [13 alone. 



It is important to note that at this order, the coun- 
terterm structure of fx/ ' f-n in a MA theory is identical 
to the counterterm structure of fx I f-n in xPT, as can be 
verified by examining Eqs. (|C1[) and (|C2I) . 



Jk 



'■K 



-ml) 



MA 



P 



MM) • 



(29) 



This can be understood with the arguments presented 
in Sec. IIIB| and the knowledge that the lattice spacing 
artifacts are flavor-blind. 



B. K K 1 = 1 scattering length, (Xkk 

Next, we discuss the form of the 1=1 KK scatter- 
ing length, for a MA theory with arbitrary sea quarks, for 
which the full functional form is provided in Appendix [El 
The two Kaon system is theoretically ideal for testing 
the convergence of SU (3) xPT, however experimentally 
much more difficult to study. But recent progress with 
lattice QCD simulations has allowed the 1=1 KK sys- 
tem to be explored within the MA framework. Thus one 
can use lattice QCD in combination with the appropri- 
ate MA effective field theory to explore the convergence 
of SU(3) xPT [2l|, or whether a generalized version of 
xPT is a more appropriate description of nature [70] . In 
fact it has only recently been confirmed that the stan- 
dard SU{2) yPT power counting is phenomenologically 
correct jsjj EH, H2] , by comparing our theoretical knowl- 
edge of the two-loop tttt scattering [83j, [84( , the pion scalar 
form-factor (85| and the Roy equation analysis (86[ with 
the recent experimental determination of the pion scat- 
tering lengths 0, [1|. 

The 1=1 KK system has several features in common 
with the / = 2 7T7r system discussed in Ref. [34|. Firstly, 



the 1=1 KK system does not have on-shell hairpins in 
the s-channel loops. Secondly, the scattering length does 
not depend upon the mixed valence-sea mesons when ex- 
pressed in terms of the lattice-physical parameters, and 
finally, the only counterterm at NLO is the physical coun- 
terterm of interest. The form of the 1=1 KK scattering 
length is given by 



' K 



'K 



'K 



' 111 1 V 



C.. In 



+ Q,- 32(4 7 r) 2 L^( A1 ) 



(30) 



where the various coefficients, C^, are provided in 
Eqs. (lE3l)-(lE6l). 

One important point is that the counterterm for the 
1 = 1 scattering length, L 1 ^ is identical to the I = 2 tttt 
scattering length counterterm, 



T i=i 



rl=2 

7T7T 



2Li + 2L 2 + L 3 - 2L A - L 5 + 2L 6 + L 8 . (31) 



Before discussing this scattering length in more detail, we 
first give the result in xPT> as this has not been presented 
in the literature to the authors' knowledge. 




with given in Eq. (|31[) , and we have used the leading 
order meson mass relations to simplify the form of this 
expression. 

The equality of the I = 2 tttt and 1 = 1 KK scat- 
tering length counterterms allows us to make a predic- 
tion for the numerical values of txiko^kk one should ob- 
tain in a simulation of this system with domain- wall va- 
lence quarks on the MILC configurations. To do this, we 
must first convert the counterterm, obtained by 

NPLQCD in Ref. [lj| from the effective theory with two 
sea flavors to the theory with three sea flavors. For PQ 
and MA theories, there is an additional subtlety which 
arises in this matching. If we match the xPT forms of 
m^a 1 ^ 2 in SU(2) to SU(3), then we arrive at the equality 
(with the conventions defined in Appendix [Dj) 




(33) 
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This leads to an exact matching between the SU(2) and 
5/7(3) theories, in which all of the strange quark mass 
dependence at this order, which is purely logarithmic, is 
absorbed in the 5/7(2) Gasser-Leutwyler coefficients |38| . 
If we naively attempt to match the 5/7(4|2) to 5/7(6|3) 
MA/PQ expressions for m^a 1 ^ 2 , using Eqs. (|D1|) and 
(|D2|) . one arrives at the relation 

^ 2 (M) = 32(4 7 r) 2 L^( A1 )-iln^)-i 

-E -T Fn(ml/m 2 x ), (34) 

where the functions, T n {y) were first determined in 
Ref. [34| and are given in Eqs. (|D3[) . All of the new 
terms in this matching arise from the extra hairpin in- 
teractions present in the (6 13) theory which are not 
present in 5/7(4|2). One can show that these terms are 
formally higher order in the 5/7(4|2) chiral expansion, 
but nevertheless we will see that they are not negligible. 

The NPLQCD collaboration has recently computed 
m^a 1 ^ 2 and used the 5/7(2) extrapolation formula to 
determine l 1 ^ 2 19]. Adjusting for conventions and in- 
cluding their largest uncertainty, they determined 

^(4^)- -10.9 ±1.8. (35) 

Starting with this determination, we can then compare 
the hairpin contributions in 5/7 (4 1 2) to those of 5/7 (6 1 3) 
and also to the physical contribution at NLO. We collect 
these results in Table [H 

The two- flavor hairpin effects, listed in row (b) of Ta- 
ble HI are not small relative to the (scale independent) 
xPT NLO contributions (a), and for the lightest two 
masses shown, are of the same order of magnitude. How- 
ever, when we consider the three flavor theory, we see 
that the additional hairpin effects, (c), are of the same 
order as the two flavor hairpin contributions (which also 
contribute in the three flavor theory), but opposite in 
sign. Taking into account all of the hairpin contributions 
by using the three flavor theory, one observes that the 
sum of these unphysical effects, (d), is approximately an 
order of magnitude smaller than the physical NLO effects 
of xPT, but increases in importance as the pion mass is 
reduced. This justifies our assumption of the determi- 
nation of l 1 ^ 2 given in Eq. ([531) . This also justifies the 
SU(3) — » 5/7 (2) matching given in Eq. ([33]) . and explains 
the success found in Ref. [lj| of using the 5/7(2) yPT 
formula to determine l{^ 2 , as the unphysical hairpin cor- 
rections to this formula provide a relative shift of about 
10% to the NLO contributions for the masses simulated, 
which is roughly the size of their largest quoted error. 

In Fig. [21 we plot the absolute values of the various 
NLO contributions to the / = 2 irir scattering length as 
a function of the pion mass, which highlight the impor- 
tance of these hairpin effects. Their relative importance 
is enhanced for the values listed in Table [J because of the 
large cancellation of the counterterm and chiral log at 



1 = 2 7T7T PQ effects for 2 and 3 flavors 

0.02 I 1 1 1 1 ! 1 1 1 1 ! 1 1 1 1 1 1 1 1 1 1 1 




100 200 300 400 500 

m„ (MeV) 



FIG. 2: We plot the absolute values of the various NLO con- 
tributions to m^a!^ 2 listed in TableQ] The NLO xPT contri- 
bution is given by (a) (green), which demonstrates the large 
cancellation of the counterterm and the chiral log for light 
to medium pion masses. The long-dashed curve (red) is the 
2-sea flavor hairpin effects, (6), which are the same order of 
magnitude as (a), for < 400 MeV. When the new 3-sea 
flavor hairpin effects, (c) (blue), are added to the 2-sea flavor 
effects, one finds that the total 3-sea flavor hairpin effects, (d) 
(black), are small compared to (a) for > 250 MeV. 



NLO, (a). It is clear that these effects will become more 
important as one moves further into the chiral regime, 

(m„ -> 0). 

Given the small contribution of the NLO hairpin effects 
to m T a^ 2 , we can use the determination of l^ 2 (fi) in 
Eq. ([35]) [l9l | , and the matching of Eq. ([33l) to determine 
and thus predict values of mjj-a^, which we 
provide in Table [TTJ We provide both the comparison of 
the NLO effects as predicted by both the MA theory as 
well as 5/7(3) xPT, which we compare to the tree level 
prediction, as well as the total scattering length through 
NLO. We find that similar to m^a^ 2 , the NLO hairpin 
effects for mxa 1 ^ are only about 10% of the NLO ;\PT 
value, less than the accuracy we claim here. We find that 
a current MA lattice determination of mKa 1 ^^ will not 
be sensitive to the unphysical hairpin contributions with 
the expected level of uncertainty, as can be seen by the 
predicted MA and SU(3) values. However, for both the 
MA and 5/7(3) theories, the NLO contributions are 15- 
30% correction to the LO term showing a convergence 
expected by power counting. 

The first error is due to the uncertainty in the determi- 
nation of L 1 ^- we obtain from the matching in Eq. ([34|) 
and the extraction of Z^= 2 from Ref. [19], Eq. (J35J). This 
uncertainty includes estimations of the two-loop contri- 
butions to m^a 1 ^ 2 in SU(2) xPT. The second uncer- 
tainty listed in Table [TT] is a power counting estimation 
of the NNLO contributions to rriKa 1 ^. Some of these 
effects are already included in the first uncertainty but 
a conservative estimate of our predicted error is to add 
these uncertainties in quadrature. 
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TABLE I: Hairpin contributions to m n a^ 2 . We provide the various hairpin contributions to the 7 = 2 -ktt scattering length 
for both the 2-sea flavor, (6), and 3-sea flavor theory, (d), which we compare to the xPT NLO contribution, (a) and LO 
contribution, top-row. In row (c), we give the new hairpin effects which arise in the 3-flavor theory, and in (d) we provide the 
total 3-sea flavor hairpin effects. 



(a) 

(b) 
(c) 
(d) 



m n (MeV) 




2-irm* 

'(47TJV)' 



[ 3 ln(^)-l-^)] 



WTT) 4 



(i n> i 



2irm% 



_2?rmJ 
W77) 1 



E 4 „=i(%)"^(^/m3 

>ai+Et 1 (if)"^(^/^)i 



293 



354 



493 



592 



-0.156 -0.218 -0.372 -0.483 

0.00460 0.00140 -0.0314 -0.0818 

0.00359 0.00327 0.00254 0.00207 

-0.00243 -0.00289 -0.00371 -0.00396 

0.00116 0.00040 -0.00117 -0.00188 



TABLE II: Predictions of mKa^K- We use the equality of the rriKa^KK an d m n a^ 2 counterterms (expressed in lattice-physical 
parameters) to predict the values of rnxOxx which would be computed in a MA lattice formulation with domain wall valence 
quarks on the MILC staggered sea quarks, which we compare both to the tree level prediction as well as the SU(3) 
prediction, for values of ttik/ fx taken from Refs. [l9l.[23|. We also provide a prediction of the scattering length at the physical 
point. The first error is due to the uncertainty in the determination of Lfof from Eq. (|34[) and the value of li^ 2 , Eq. (|35[) 
determined in Ref. [Tj|. The second error is a power counting estimate of the NNLO contributions to the scattering length. 



m K : f K (MeV) 


577 : 172 593 : 171 639 : 173 690 : 177 


mKa 1 ^ (NLO: MA) 
m K a^(NLO: SU(3)) 


-0.447 -0.479 -0.542 -0.605 

-0.091 -0.113 -0.162 -0.223 
-0.084 -0.107 -0.157 -0.217 


mKogx(MA) 
m K a I £^{SU{3)) 


-0.540 ± 0.069 ± 0.026 -0.592 ± 0.079 ± 0.031 -0.704 ± 0.102 ± 0.048 -0.828 ± 0.127 ± 0.072 
-0.531 ± 0.069 ± 0.026 -0.586 ± 0.079 ± 0.031 -0.699 ± 0.102 ± 0.048 -0.823 ± 0.127 ± 0.072 


physical point 


496 : 161 


m K aM{SU(3)) 


-0.424 ± 0.049 ± 0.012 



C. Ktt I — 3/2 scattering length, a K7T 

The Ktt system is also an interesting laboratory for 
exploring the three flavor structure of low energy hadron 
interactions, and moreover it is experimentally acces- 
sible with proposed studies by the DIRAC collabora- 
tion [89j. There has recently been a direct MA lattice 
QCD determination of the I = 3/2 Ktt scattering length, 
which in combination with the theoretical knowledge of 
the NLOyPT I = 1/2 and I = 3/2 Ktt scattering 
lengths [SO, [HI, [H, [93| has allowed a determination of 
both isospin scattering lengths |25j |. There is addition- 
ally a two-loop computation of Ktt scattering in SU{3) 
xPT which studies the convergence of the theory with 
standard power counting [94| . Before embarking on a 
study of the two- loop effects with lattice QCD, one must 
first understand the lattice corrections at NLO. This is 
the motivation for this section. 

The tree level 7 = 3/2 Ktt scattering length is given 

by 



/ . n 7=3/2 

[m 7r + m K )a K ^ 1 =- 



1=3/2 



m^rriK 



or 



(36) 



where fiKir is the reduced mass of the Ktt system. We 



chose to express our extrapolation formulae in terms of 
the product fx fir since this symmetric treatment of the 
K and tt mesons provides the simplest form of the scat- 
tering length. We find, however, that the I = | scatter- 
ing length still depends on the mixed valence-sea meson 
masses, and therefore on the parameter Cmix- Conse- 
quently, accurate chiral extrapolations of this scattering 
length will require a determination of the value of Cmix 
appropriate to the particular sea quark discretization 
used in the simulation. The form of the MA 1 — 3/2 Ktt 
scattering length is 



7=3/2 



+ 



f l Kir 
^f K fir 

"(m K - 



32m K m % I=2 

1 ■ ■ IT IT W 



fK fir 



L5M 



Kir, 3/2/ 



(AO , (37) 



where am?' 3 (jjl) is the valence- valence (including 
valence-ghost) contribution to the scattering length and 

a^ 7r,3 ^ 2 (/i) is a non- vanishing contribution from mixed 
valence-sea mesons to the scattering length. The other 
important thing to note is that there are two countert- 
erms for this scattering length which can both be de- 
termined through its chiral extrapolation formula, but 
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can also independently be determined in other processes; 
L 5 (/i) can be determined independently by /k/Itt an d 
L-mt {/J*) can be determined either with 1 = 2 nir or 
1 = 1 KK scattering. Now we have explicitly demon- 
strated that the four observable quantities, a^ 2 , a 1 ^, 

a K~^ 2 an( i fi</ fm when expressed in terms of the lattice- 



physical parameters, only share two linearly independent 
counterterms through NLO in MA (and PQ) %PT. 

Before continuing, we provide the continuum SU (3) 

1—3/2 

xPT form of , which is the same as can be con- 

structed from Ref. |93| with the NLO shift of fx — > f*, 



1=3/2 
(J>Kir a K Ti - 



+ Ktan arctan 



K n In ( --j- I + K K In I "7^2 J + Kl ' 111 I "TTi" I ~~ "g* TO A"V 



2(to^ - m 7T )y / 2m 2 K + m K m^ - m 2 \ 
(2mR- - TO 7r )(TOK + 277V ) i 

32771^777^ r J =2 / s 

7 7 7T7T W 

/if Jtt 

I 



Ik fx 



L^)\, (38) 



with 



m 2 \\m 2 K + 22ttixTO7t — 



(39) 



mif 134m| 



9 m K + 55771^771^ — 16m^. 
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A" 



-SGm^ — Ylrrv^m^ + irixm 2 + 9m^ 



16toa"TO w ^2m 2 K + niKm^ — m 2 



m,K — fn-K 



(40) 
(41) 

(42) 



In Appendix [0 we provide the full form of the MA 
7 = 3/2 Ktt scattering length. Here we wish to examine 
the valence-sea contribution in more detail. This contri- 
bution to the scattering length is given by 20 



C Fs In 



E 



l Fs 



C F d In 



F=j,l,r 

2 > 
Fd 



4mKm 7r J(m 2 Fd ) 



(43) 



20 We note that the summation over sea flavor, F, implicitly in- 
cludes the appropriate factors for staggered sea quarks, the sum 
over taste and the factors of N s = 1/4 which arise from the 4 th - 
rooting trick. For other brands of quarks, this is simply a sum 
over the sea flavors, j, I and r. 



where 

Cfs - 

CFd - 

J(M) : 



Amr K m v - m 2 Fs (m K + m OT ) 
ArriKm 2 - m 2 Fd (m K + m n ) 



m,K — 



(44) 
(45) 



,VM 2 



arctan 



(niK - m^y/M 2 - to 2 \ 
M 2 + rriKm^ — m 2 J 

(46) 



The scale dependence in a^ s n ' 3 (fi) can be shown to be 



n 



Air, 3/2 



E 

F=j,l,r 



\n(n 2 )(m K ~ m w f 



(47) 



and as claimed, independent of both the lattice spacing, 
a, and the sea quark masses, and is absorbed by L$(fi). 
We stress again that the 7 = 3/2 Kir scattering length 
depends upon mixed valence-sea mesons, which receive 
lattice spacing dependent mass shifts proportional to the 
unknown quantity, Cmix- This quantity is currently un- 
known for all variants of MA lattice QCD and must be de- 
termined for a correct extrapolation of MA lattice QCD 
simulations. For this reason, we do not provide a table 

1=3 12 

with post-dictions of HK^a Kv 1 . 



IV. DISCUSSION 

Mixed action simulations provide a promising solution 
to the problem of performing fully dynamical simulations 
with light quarks which are under theoretical control. Re- 
cently, several simulations have been performed using the 
publicly available MILC staggered lattices with domain- 
wall valence quarks and the future for simulations us- 
ing Wilson sea quarks and GW valence quarks seems to 
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be bright. This work shows that mixed action simula- 
tions with Ginsparg- Wilson valence quarks are theoret- 
ically clean. We have shown that the counterterms ap- 
pearing in mesonic scattering lengths are precisely those 
that occur in QCD, so that one can in principle mea- 
sure these counterterms with a single lattice spacing. We 
also find that the same chiral extrapolation formulae can 
be used to describe mixed action simulations with GW 
quarks with mild restrictions on the type of sea quark 
discretization used — provided, of course, that QCD is re- 
covered in the continuum limit. Thus, our results hold for 
simulations with domain wall and overlap quarks, Wilson 
quarks (0(a) improved and twisted mass quarks at max- 
imal twist) as well as simulations using rooted staggered 
quarks (assuming that the 4 t/l -rooting procedure is valid 
and that the replica method correctly captures all of the 
non-locality introduced by the rooting procedure - which 
has been argued to all orders in perturbation theory [95|). 

It was previously observed in Ref . [34| that the irir scat- 
tering length does not depend on the parameter Cmix of 
mixed action chiral perturbation theory. Here, we find 
that this also holds for the KK scattering length. How- 
ever, the Ktt scattering length does depend on Cmix, and, 
therefore, accurate chiral extrapolations of mixed action 
data will require a measurement of this quantity. How- 
ever, we have also computed the ratio Jk/ f-n in mixed 
action chiral perturbation theory, which depends upon 
on Cmlx- By varying Cmix over a broad range of values, 
we find the impact to be modest, on the order of 5%. 
In addition, taking into account the small hairpin cor- 
rections to m^a^ 2 [HI and the equally small predicted 
corrections to ttiko^kk^ Table HT1 we expect the impact 

1—3 /2 

of Cyiix on a Kn to also be small at this order of preci- 
sion, 0{e^ n ,e 2 n e 2 a ,e\). In Table HT1 we provide predictions 
of mK^KK f° r various values of rriK and also provide a 
prediction at the physical point. 

In Sec. Ill BL we have demonstrated why the use of 
lattice-physical parameters (or on-shell renormalization) 
significantly simplifies the form of the extrapolation for- 
mulae for mesonic systems. We stress that these argu- 
ments do not depend upon the momentum of the system, 
nor upon having only two external mesons, and thus will 
be applicable not just for scattering lengths, but also for 
other scattering parameters, such as the effective range, 
as well as for N > 2 mesonic systems. In the appen- 
dices we have provided explicit NLO extrapolation for- 
mulae for the meson masses and decay constants as well 
as the three scattering lengths discussed in this paper, 
for arbitrary sea quark discretization schemes, expressed 
in terms of the PQ parameters we introduced in Eq. (fT4|) . 
A thorough understanding of the lattice spacing effects 
at this order will require knowledge of the counterterms 
in the masses and decay constants. 

We would like to conclude with a small point and a few 
suggestions. If one is interested in removing the unitarity 



violating effects in MA lattice simulations, for the low- 
energy dynamics of the system, then theoretical analysis 
unambiguously advocates the tuning A rs = Aj U = 0, 
which is the generalization of m qsi!a — m Qval — for PQ 
theories. This is the most QCD-like scenario for MA the- 
ories in which the unitarity violating double pole prop- 
agators in Eq. (TT2|) are tuned to zero. It has recently 
been shown that this double pole structure of the flavor 
neutral propagators persists to all orders in PQxPT [96| . 
and thus this will be the appropriate tuning to higher 
orders as well. From the point of view of doing chiral 
physics, this is not desirable for the coarse MILC lattices, 
as the lattice spacing shift to the taste-identity staggered 
mesons is a 2 Aj ~ (450 MeV) 2 , which would make for 
heavy pions. Therefore we caution users of MA lattice 
simulations to remember the existence of these unitarity 
violating effects present in current MA simulations. 

The simplified form of MA/PQ extrapolation formu- 
lae for the two-meson systems is particularly dependent 
upon the implications of the chiral symmetry of the va- 
lence quarks. However, we conjecture that a similar, but 
not as strong, simplification will occur for other hadronic 
observables, in particular for nuclear physics as well as 
heavy meson observables, if also expressed in terms of 
lattice-physical parameters, which will lead to improved 
chiral extrapolations. This is supported b y th e recent fits 
of the NPLQCD collaboration [H HUTU, HI H| and 
the LHP collaboration [97]. Based upon our theoretical 
understanding of effective field theories designed to in- 
corporate lattice spacing artifacts, we expect that even 
for fermion discretization schemes which do not have chi- 
ral symmetry, the use of lattice-physical parameters (on- 
shell renormalization) will in general simplify the chiral 
extrapolation formulae and improve chiral fits. 
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APPENDIX A: m„ AND f v FOR 2 SEA FLAVORS 

In this Appendix, we provide the explicit formulae for 
the pion mass and decay constant in a two-sea flavor MA 
theory. These were first computed in Refs. [HI, [3(|. Here 
we provide the answers expressed in terms of the PQ 
parameters we introduced in Eq. (fl"4"|) . 
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APPENDIX B: MESON MASSES 



In this Appendix we collect the pion and kaon mass and decay constant for a three-sea flavor MA theory. These 
were first computed in Refs. (l5l . [30| . Here we provide the answers expressed in terms of the PQ parameters we 
introduced in Eq. (fl4|>. 
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Note that the lattice spacing dependent counterterms for the meson masses have the same coefficient. This is because 
the discretization scheme is flavor-blind. 



APPENDIX C: DECAY CONSTANTS AND f K /.U 
The pion decay constant is given by 
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while the kaon decay constant is 
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The two important things to note are that the addi- 
tive lattice spacing modifications to the decay constants 
are the same and also that at this order, they can be ab- 
sorbed into a redefinition of the Lagrangian parameter, 
/. We can then use these formulae to estimate the size 
of the corrections to the recent determination of L 5 (pi) 
by NPLQCD 23]. Thus, we form the ratio 



where, using Eqs. (|Cip and (|C2[) . and the tuning used 
in Ref. [23j which was to set the valence-valence meson 
masses equal to the taste-£s sea-sea mesons, we have 
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APPENDIX D: n + n+ SCATTERING 



Three sea Quark Flavors, m,a^ 2 



For completeness, we provide the formulae for the I = 
2 7T7T scattering length determined in both MAvPT and 
PQxPT in Ref. [34], for both two and three flavors of sea 
quark. 



Two sea Quark Flavors, m^a 
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APPENDIX E: K+K+ SCATTERING 



where 



The I — 1 KK and 7 = 3/2 Kir scattering lengths 
involve lengthy expressions. Therefore, we introduce the 
following notation to make the answers more presentable. 



m K = km^ , Aj u = S^m^ , A rs = S r 



The I ~ 1 KK scattering length is given by Eq. (|30|) . 
which we repeat here for convenience 
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APPENDIX F: K + n + SCATTERING 



The Kit scattering length at 7 = 3/2 is given by: 
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where a^ T,3 ^ 2 (fi) is given in Eq. ((43 
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and the coefficients, Cp<j s and the function, J{m) are defined in Eqs. (pT4"|) - (f4"6")) . We reiterate that the ln(/i 2 ) 

dependence in a^ 7r ' 3 ^ 2 (/i) only depends upon the valence- valence meson masses, Eq. (|47[) . as we argued in Sec. Ill B] 
The valence-valence (and valence-ghost) contribution to the scattering length is given by 
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We use the notation defined in Eq. (|E1[) to simplify the form of these coefficients. We find 
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9fc 6 - 126/c 5 + 113fc 4 + fc 3 (139 - 100(5 2 J + fc 2 (-167 + 64(5 2 J + fc(41 + 36<5 2 s - 20(5 4 s ) - 9 



+ 25* 



108/c 8 - 488fc 7 + 3fc 6 (73 + 185 2 s ) + fc 5 (570 - 828<5 2 s ) + 3fc 4 (-47 + 274<5 2 J + 81 - 54(5 



6fc 3 (36 - 975 2 s + 625 4 J + 3fc 2 (-89 - 2U5 2 S + 112(5 4 J + fc(134 + 665 2 s + 36(5 4 s - 40^ s ) 



+ 85? 



108fc 



10 



56k 9 + 2fc 8 (-251 + 54(5 2 J - 4fc 7 (-87 + 865 2 s ) + 3fc 6 (232 - 235 2 s + 9(5 4 J 



- 6fc 5 (82 - 104<5 2 s + 57(5 4 s ) + 3fc 4 (-124 + 13<5 2 s - 
+ fc 2 (124 - 159(5 2 s - 141(5 4 s + 88<5^) + fc(36 + 80(5 2 



fc 3 (164 - 360(5 2 



273(5; 



1125 r 6 s ) 



57(5 4 s + 24(5 r 6 s - 105 8 s ) - 27(2 - 3<5 2 s + 5 A rs ) 



+ 8{2k 2 -2 + 5 2 rs ) 2 36fc 8 - 24fc 7 + fc 6 (-85 + 18<5 2 J + fc 5 (40 - 285 2 rs ) + fc 4 ( 71 - 46(5 2 J + 9 
- 2fc 3 (4 - 33(5 r 2 s + 8(5 4 J + fc 2 (-31 + 10<5 2 s - 4«5 4 J - 2fc(4 + 5 2 rs 10«5 4 s + 25 e rs ) ' - - '- 



(F5) 
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A = 



: {4096fc(fc 2 - l) 6 + U8 2 rs {k 2 - if {9k 2 + 56fc + 9) 



9(fc 2 - 1) 2 (2 - 2fc 2 + S] u + 2<5 2 s )(4fc 2 - 4 + S 2 U + 25 2 s f 
+ 96S 4 s (k 2 - l) 4 (9fc 2 - 8fc + 9) + 16(fc 2 - l) 3 (27fc 2 - 88k + 27)Sf s + 8{k 2 - if (9k 2 - 40k + 9)S* ( 



3 " 



32(fc 2 - l) 5 (9fc 2 + 14k + 9) + 96S 2 s {k 2 - if {3k 2 - 35k + 3) - 24<5 4 s (fc 2 - l) 3 (9fc 2 + 166fc + 9) 



- 8<5^(fc 2 - l) 2 (36fc 2 + 155fc + 36) - 726* s (k 4 + fc 3 - k - 1) 



48fc(fc 2 -l) 4 (6fc 2 + 5) 



+ 12S 2 Jk 2 - l) 3 (12fc 3 - 27fc 2 - 136fc - 27) - 6<5 4 s (fc 2 - l) 2 (36fc 3 + 63fc 2 + 230fc + 63) 



- 36<5^(5fc 5 + 3fc 4 + 2fc 3 - 7k - 3) - 36<5^fc 3 



2(fc 2 - l) 3 (36fc 3 + 27fc 2 + 394/c + 27) 



+ 2S 2 s (k' 2 - l) 2 (108fc 3 + 63fc 2 + 323k + 63) + 185* s (+9k 5 + 3k 4 - 4k' A - 5k - 3) + 365^ s k 3 



(fc 2 - l) 2 (36fc 3 + 9fc 2 - 22fc + 9) + 9(4fc 5 + fc 4 - 5fc 3 + fc - 1)<5 2 S + 9k 3 S 4 s } , 



(F6) 



4fc 



Han 



(fc - 1)(2 - 2fc 2 + 6% + 25 2 rs ) 2 ^J{k-lf{k + l) 



: arctan 



v / (fc-i) 3 (fc + i) 



fc 2 + fc 



f^){8(fc-l) 4 (fc + l) 3 +45 2 s (fc 2 -l) 2 



45 r 4 s (fc 3 - 2fc 2 - fc + 2) - 5 2 u [8(fc - l) 3 (fc + l) 2 + 45 2 s (fc 2 - 1) - 2S 4 S ] + S 4 U [2(k - l) 2 (fc + 1) + S 2 rs ] } 



8fc(2-2fc 2 + ^-<5 2 ^ 



- 12k + 4 - ^ - 2<52J( 4 fc2 - 4 + <5 2 „ + 26 2 s ) 



9(fc-l) 2 (2-2fc 2 + <5 2 „ + 2J 2 s ) 2 



x arctan 



/ ^(8fc 2 - 12k + 4 S 2 U - 20tf - 4 + S 2 U + 25 2 rs ) x 
4fc 2 + 6fc - 4 + W u + 2S 2 S ) 



(F7) 
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